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INTRODUCTION 



The general validity of the superposition principle of quantum mechanics is questioned by an 

J], and is subject to an intense experimental verification 5- 10]. 



increasing number of scientists [1 
The possibility that the su p er p osition principle might have only a limited range of validity is fore- 



seen by collapse models 



11 



221 ] . which predict small deviations from standard quantum mechanics, 



for all those cases where quantum linearity plays a fundamental role. Neutrino oscillations are one 
such a case, and the goal of this article is to present the theoretical analysis and estimate of the 
effect of spontaneous collapses on the oscillatory behavior of neutrinos. 

n 

A previous analysis of this kind was proposed in [23J, based on the Penrose model of gravity 



induced collapse 



25j. This model however lacks a fu 



previous attempts to fill this gap 



ly consistent dynamical 



261 ] have been criticized 28[]. Moreover, as shown in 28(, the model 



equation, and 



fails when applied to single constituents (protons, electrons, ...), since in this case its predictions are 
in conflict with known experimental data. Therefore, the application of gravity induced collapse 
models to neutrino oscillations is rather delicate. 

Here we will com put e the spontaneous collapse effect on neutrino oscillations using the mass 
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proportional version 
together with the GRW model 
equation: 



121 of the CLS model 



Q. 



131 ]. which is widely used in physical applications, 



Its dynamics is described by the following stochastic differential 



d\4>t 



- % -Hdt + / dx (M(x) - (M(x)» dWtU) 
n too 



7 
2m\ 



dx (M(x) - (M(x))) 2 dt 



(1) 



where the operator H is the standard quantum Hamiltonian of the system and the other two 
terms induce the collapse of the wave function in space. The mass too is a reference mass, which 
is taken equal to that of a nucleon. The parameter 7 is a positive coupling constant which sets 
the strength of the collapse process and (M(x)) = {fy |M(x)| cj>t), where M(x) is a smeared mass 
density operator: 



M ( x ) = ^Z m i^Z j dy9 ( x ~ y ) ^] ( y > s ) ^3 (y> s ) 



(2) 



if)j (y, s), if)j (y, s) being, respectively, the creator and annihilation operators of a particle of type 
j, having mass to, and spin s, in the space point y. The smearing function g(x) is taken equal to 

, . 1 



9W 



-x 2 /2r2. 



2irrc)' 



(3) 
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where rc is the other new phenomenological constant of the model. Finally, Wj (x) is an ensemble 
of independent standard Wiener processes, one for each point in space. The standard numerical 
value of the correlation length rc is [3]: 



rc ^ 10 5 cm, (4) 



rave been proposed. The 



3: 



while, in the literature, two different values for the collapse strength 7 
first value has been originally proposed by Ghirardi, Pearle and Rimini 

7 ~ l(T 30 cm 3 s _1 (5) 

in analogy with 1 lie GRW model [14f |. The second value has been proposed by Adler, inspired by 
the analysis of the process of latent image formation according to collapse models, and amounts 
to Q: 

7 ~ l(T 22 cm 3 s _1 (6) 



Working with non-linear equations is notoriously difficult. As shown e.g. in [29], the experimen- 
tally testable predictions of the model — when averaged over the noise — do not change if the real 
noise Wt (x) is replaced by an imaginary noise iWt (x). In this way, one loses the collapse proper- 
ties of the equation. However, the advantage of having an imaginary noise is that the evolution is 
described by a standard linear Schrodinger equation with a random Hamiltonian: 

H TOT = H - VtX^E / dyw(y,t)^(y, 3)^,3), (7) 
j s 

where 

w(y,t) = J tfc0(x-y)6(x), (8) 

and £t( x ) = dWt{x)/dt is a white noise field, with mean equal to zero and correlation function 
E[£t(x)£ a (y)] = 6(t — s)<5(x — y). As such, w(x,t) is a Gaussian noises field, with zero mean and 
correlation function: 

EKx.tMy,*)] = 5(t- S )F(x-y), F(x) = Jl e'^^C (9) 

{yA-nrc) 

We are interested in the relativistic generalization of the Hamiltonian ([7]). The most natural 
choice is: H TOT = H D + N(t) where (in the case of just one type of particle): 

H D = J dx.Ho (x) = Jd^ (x) [-ihccx ■ V + mc 2 f3] ip (x) (10) 



4 



is the standard Dirac Hamiltonian. Here we have introduced the four- vector notation x = (ct,x), 
c is the speed of light, m the mass of the particle associated to this field, a = (a 1 , a 2 , a 3 ) with 
a l = 7°7* and P = 7 where the 7^ are the Dirac matrices, which we take in their standard 
representation: 

( 1 \ ( a 1 \ 

7°= , y= with i = 1,2,3. (11) 

\0-lJ W / 

where 1 is the identity matrix in two dimensions and a 1 are the Pauli matrices. The noise term 
instead is given by: 

N(t) = JdxAT(x) = -h y / ! ujdxw(x)^(x)^(x), 7ro = 7 Q^ (12) 



It can be showdjthat the Hamiltonian H TOT = H u +N(t) defined in this way, in the non relativistic 
limit, reduces to the Hamiltonian of Eq. ([7]). 

Working with plane waves gives rise to unphysical infinities, since they are not normalizable. 
To avoid potential problems, we use the box normalization, i.e. we confine our fields in a box of 
length L, and we impose periodic boundary conditions: ip (t,x) = tp (t, x + L), where L is a vector 
with all the components equal to L. In turn, the momentum is discretized: p = 2 -^k with k % e Z 
and z = 1,2,3 labels the spatial components. Then the Dirac field, in the interaction picture where 
we choose H D as the unperturbed Hamiltonian and N (t) as the perturbation (the noise coupling 
x/TV^ is very small), takes the usual expression^: 



2 +00 -. / .1 
1 / mc z 



^) = E E ^W^[fc(p,-)«(p,-)e-*^ w ^ + crt(p ia )i;(p,*)e*^¥ k -' e J, (13) 

s=l k=— 00 • p 



where u and v are the usual Dirac spinors, E p = y p 2 c 2 + m 2 c 4 is the energy and b and d are 
operators satisfying the standard anti-commutation relations. We also recall the relation between 
the evolution operator U(t) in the Schrodinger picture and Ui (t) in the interaction picture [30! ]: 

U(t) = e-i H ° t U I (t), (14) 

and we set the initial time to 0. 

One can question why we use a relativistic Hamiltonian in a model which is not relativistic, 
since the correlation function of the noise is not Lorentz-invariant. Our approach to this issue, 
is that collapse models are phenomenological models emerging from a pre-quantum theory yet to 



1 A proof of this can be found in section 2.3 of [341 ]. There, the proof is worked out for the Dirac equation coupled 

with an electromagnetic field. This is the same casei as ours, if one sets A (x) = and eAo (x) — -fi^/7 m w (x). 

2 Here and in the following we use the notation of [341 ]. 
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be discovered. The noise field is a real cosmological field (whose nature is yet to be investigated) 
which naturally defines a privileged frame, most likely corresponding to the co-moving frame of the 



universe. Hence we see no contradiction in analyzing relativistic phenomena with the CS 
For attempts towards a fully relativistic formulation of collapse models, one can refer to 



model. 



31 



32]. 



II. THE TRANSITION AMPLITUDE 



In the more general approach to the problem of neutrino oscillations, we consider n flavour 
eigenstates, which will be labeled by greek subscripts \ va) and each of them is a linear combination 

of n mass eigenstates that will be labeled by latin subscripts 

n 

\v<x) = y^Uqj Wj) ■ (15) 
3=1 

Here, U is the nxn mixing matrix, which relates the two different bases; since the flavour eigenstates 
are supposed to be orthonormal, as well as the mass eigenstates, U must be unitary. 

We take a neutrino in an initial flavour eigenstate, and compute the probability of finding 
it in another flavour eigenstate, after some time t, assuming that the dynamics is governed by 
the Hamiltonian H TOT . We assume that the neutrino has definite initial and final momenta. 
Mathematically, we will compute the following quantity: 

n 

T a ^p = (vp;Pf,Sf \U (t)\ u a ;pi,Si) = ^2 U "i U /3i (^;P/> S / P (t)\ (16) 

here U (t) is the time-evolution operator of Eq. (|14|) while \v a ;Pi,Si) is the flavour eigenstate a 
with initial momentum pj and spin Sj and |z/g;p/,Sf) is the final state with momentum pf and 
spin sf. Since the Hamiltonian is the sum of Hamiltonians associated to different mass eigenstates 
(H = ]Cj=i ^ i s convenient to expand the flavour eigenstates into the mass eigenstates, as we 
did in Eq. ()16|) . The form of the Hamiltonian also implies that U(t) factorizes: U (t) = YYk=i Uk (t)- 
Here Uk (t) is the time evolution operator related to the Fock space of the neutrino having a definite 
mass mfc. This is an important property, because it implies that if i ^ j: 

(vf,pf } s f \U(t)\vj;pi, si) = (ill \Ui (t) | Hi) ... (ui;pf,Sf\Ui (t)\ Oj) ... (Qj \Uj (t)\ u^p^si) ... 

...(n n \u n (t)\n n ) = o (17) 

since (i/j; pf, Sf \Ui (t)\ Qi) = 0, as one can check with a direct calculation. Therefore we can write: 

n 

T a _>p = Y, V «i V h [< n i l^i - (^;P/' S / \ U 3 vf,Ph s i) ■- ( n n Pn (*)| n„}] , (18) 

3=1 
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which reduces the entire calculation to a 1-particle computation. In the next section we will 
focus our attention on the matrix element (i/j;pf,Sf \Uj (t) \ iSj;pi,Si) since, as we will show, the 
remaining terms contribute with an unimportant global phase factor. 



III. THE MATRIX ELEMENTS 



We now focus on the main part of this work. What we need to compute the 1-particle matrix 
element: 



T(jp f ,sf,pi,sf,t) = (p f ,s f \U(t)\pi,Si) = e ** ft (pf,8f\Ui(t)\pi,Si). (19) 

Since this part of the computation is the same for every mass eigenstate, we have dropped the 
label j. We expand the evolution operator by means of the Dyson series up to the second order: 



Ui (t, 0) ~ 1 - l - J dh : N T (h) : J dti J 



dt 2 : Nj (h) :: Nj (t 2 ) :, (20) 



where Nt (t) is the interaction picture representation of Eq. (|12p and : ... : denotes the normal 
ordering^. Accordingly, the transition probability is the sum of three terms: 

T (P/, sf, Pi, sf, t) = e~i E ^ T (0) (pf, s f ; pi, sf, t) + T (1) (p f , s f ; p { , s { ; t) + T (2) (p/, s f ; p { , sf, t) , 

(21) 

corresponding to the first three terms of the Dyson series. We now give a pictorial representation of 
each term by means of Feynman diagram and we compute each of them. The first term corresponds 
to the free propagation: 

T^(p f ,s f ;pi,Si;t) = 1 * f 

where the solid line represent the particle. This term is trivial: 

T (0) (p f ,s f ;pi,Si;t) = (p f ,s f \pi,Si) = 5 SfSi 5 Pf , Pi . (22) 
The second term correspond to the diagram: 

T (1) (p f ,sf,Pi,Si;t) = 



1 V 

where the dotted line represents the noise field. This term is: 

r (1) (pf,sf,Pi,Si;t) ^VTm J dti Jdxiw{xi){pf,Sf\:^i{xi)^ I {x]):\pi,s i ). (23) 



3 As well know in Quantum Field Theory, the reason why we used : Ni (t) : in place of Ni (t) is that, with this 
prescription, we can remove all divergent contributions coming from tadpole diagrams. This type of divergences can 
be absorbed through a renormalization procedure, without giving any physically observable consequence [35I. 
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In order to compute the matrix element in Eq. (J23J), we use the series expansion of the fields as 
given in Eq. (|13j) . The non-null terms are those containing two b and two operators. After some 
calculations, one finds that: 

(Q\b(pf,3f) : 4>i (xi)ipi (xi) : 6 1 " (pi,Si) = 

2 +00 1 2 

= E E Tg-S^^-^u^^ 

S /=lp,p'=-oc L V^P' 

= J^— e^ p f- p ^u{ Pf ,s f ) u (ft, 8i ) . (24) 

Here we introduced the four momentum j/* = (E p /c, p). If we substitute Eq. (I24p in the definition 
of TW wc get: 

(pf,sf,Pi,Si;t) = i^/^-^=u{pf,Sf)uipi,s i ) T ^ / cfti / dxiw (x x ) e» (p/ - ^) 3 ^. (25) 

y&ihf L Jo J 

The last term in Eq. (|21|) is the more complicated to compute, since it involves the product of 

four fields. It gives the following contribution: 



T {2) (pf,8f,Pi,8i;t) = —jp(Pf,8f 



t rt 



1 



dh / dt 2 : iVj(ti) :: N z (t 2 ) : 
Jo 



Pi,Si)= (26) 

= --y y dt 1 dt 2 J dx 1 dx 2 w (xi)w (x 2 ) (p/,s/|T [: t/>j (xi) V/ (xi) :: ?/>/ (x 2 ) (^2) :] |pi, s^), 

here "T" is the time-ordering product. Using Wick's theorem, and discarding all tadpole terms, 
which involve a contraction between two fields at the same spacetime point, we have: 

T [4>la'<Pla1p2b1p2b] = ■ ^\a^\a^2b^2b ■ ~S a b (xi ~ X 2 ) S ba (x 2 - X X ) 

+iS ah (x\ - x 2 ) : tp la ip2b ■ -iSba (x 2 - xi) : ViaV^fe (27) 

where a and b label the spinor components and the Dirac propagator is: iS a b (x% — x 2 ) = 
(f2|T [V>iaV>2b] |^)- Here, Einstein's summation convention is used for the spinor indices. We 
momentarily drop the pedex / related to the interaction picture and we write the dependence on 
x\ and x 2 simply as 1 and 2 . The diagramatic representation of the different terms in Eq. (|27|) is: 




L ~*~~ 2 \^ f = 1 ^^Ia1p2b1p2b ■ , { / -.«?„/, (Xj. - ,r 2 ) >%, (.''2 - -''J ) ■ 



>— *~\. = iS ab (xi - x 2 ) : ipi a ijj 2 b : , J-^—K = ~ iS ba (x2 - xi) : tpi a ip2b ■ ■ 

We can easily see that the first term is zero, since we are studying the case with only one particle 
in the initial and final states. Regarding the second diagram, an important issue here arises. This 



diagam represents a vacuum fluctuation term, which is divergent. As well known |36(, all vacuum 
fluctuations diagrams of any order sum up to a phase factor (f2 \ Ui (t, 0)| Q), and all divergences 
cancel each other. Therefore we can write: 

(pf,8f\U I (t,0)\p i ,8i) = {Sl\U I (t,0)\n)'{pf,8f\Ui (t, 0)| Pi , s t ) cxt , (28) 

where {pf,Sf\Ui(t,0)\pi,Si) t denotes the contribution from diagrams with external fermionic 
legs. This vacuum fluctuation term is important because, together with those of Eq. (|18p . it gives a 
global phase YYk=i {^k Wki (t, 0)| Qk) independent of j, which factorizes out of the sum. Therefore, 
such terms are physically unimportant, the only relevant part being (pf,Sf\Ui(t,0)\pi,Si) t . 
From now on, we will work only with diagrams with external fermionic legs, and we drop the 
pedex "ext". 

Coming back to Eq. (|27p . we can now focus our attention on the third and the fourth term, 
that correspond to the last two diagrams. Since : ip 2b ipia '■ = — '■ ipia^b ■ for fermions, these two 
terms give the same contribution. The Dirac propagator reads: 

2 +oo 2 
s=l p=— oo P 

-9 (i 2 - h) e~kP^-^)v b 0, a) v a (p, s)} , (29) 

while the matrix element gives: 

/ — .1 mc 2 i/ T1 i-i 1 ^ . . , . ,„„. 

( P/ , S/ | : M 2b : |p i)Si ) = __— ^eftf/^e *Pi x ^u a (p/, s f ) u b ( Pi , s;) . (30) 

In the following we will need only the case p/ = Pi and st = Sj. In this case, we have: 

1 2 2 +°° 1 2 

C21 1 mc \ -> \ - 1 rac 

T ( [Pi, sf, pi, Si;t) = — 7 m — g— - 2^ 



L 3 Ei ^ ^ L 3 £L 

I s =l p=-oo p 

I I dt 2 f dx x / dx2'u;(xi)u;(x2)e^ p ^ pM )( :E ^- a;2 ^ 

JO J J 

• ^afe, Si)Ua(p's)uh(p's)u b (pi, Si) 

-J dt 2 dti J dxxj d^ 2W { Xl )w{x 2 )e^- p ^ x ^- X2 ^ 
■ u a (pi, Si)v a {p, s)v b (p, s)u b (pi, Si) I . (31) 
Using the standard relations j^ : 

\2u a (p,s)u b (p,s) = ( P > + mc j Y]v (p,s)u b (p,s) = (^-^ — —J , (32) 

jri V 2mc y afe ^ V 2mc 
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we can see that the terms containing p^j^ cancel each other, while those containing the mass give 
a 6 a b- Thus, if we also use: u(pi, Sf)u(j)i, Sj) = S Sf)Si , we obtain: 



T (2) (pi,Si;pi,Si;t) = -j m j-^— ^ 



1 mc \ - 1 mc r ' 



p=— oo 

Xe H(pf-P M )(^lM-^2 M ) > 



L 3 



p jo Jo 



dt\ I dt2 I dxi / dx2W (x\) w (#2) 



(33) 



Now we have all the elements we need, in order to compute the transition probability. We will do 
this in the next section. 

IV. THE TRANSITION PROBABILITY 



The physical quantity we are interested in, is the transition probability, which corresponds to 

2 

|Ta->0| j averaged over the noise, and integrated over the final momentum and polarization states: 

(34) 



s f Pf = ~ 00 fc=lj=l 



where: 



+00 



p k j (Pi,si]t) = ^2 eh{E<t ' ' B/t,) ' E K (p/,s/;Pi,si;*) r j (p/,s/;Pi,^;t)] , (35) 

Sf p^=— oo 

and Tj (py, sy; pj, Sj; t) is given by Eq. (I2ip . where now we have explicitly indicated the label j 
associated to the mass eigenstate rrij and E^p = y/p^c 2 + m 2 c 4 . When averaging, one has to 
remember that only terms containing an even number of noises survive (in the Feynman represen- 
tation, all products of diagrams containing an even number of dotted legs). Using this fact, and 
exploiting the Kronecher deltas of (pf,sf,pi,Si]t) (see Eq. ([22]) ). we can write: 

P kj (pi, Si ;t) = eii^-E^t + (p i)Si;t ) +lf ( Pi , Si ;t) + J< 2 >* (pi, Si ;t)] , (36) 
where we have defined: 

■AX) i + \_V^ iiE^-E^-E^p+E^ht-^ Ml)* / *\rr(X)t W 

i)k (Pii^t) = 2^ e E [ T k (p/^/sPj)^;*)^- (p/,*/;p<>*i;<)J > 

P/=— oo 

jj 2) (ft, Si ; t) = E [if } (ft, Sl ; Pi , . (37) 

We focus our attention on ij^ (pj, sf,t). Using Eq. ([25]) . keeping in mind the spinor relation 
(llfUi)* = UiUf, and performing the average over the noise, which brings in a Dirac delta in time 
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which cancels one of the two time- integrals, one obtains: 



E 



T% (pf, s f ; pi, sf, t) T^> (p f , sf, pi, Si; t 



(i) 



mjrrikC 



u(p ( p ',s f )u(pf> ,Si)u{pf> ,Si)u(p ( j> ,Sf) 



Si) 



(k) 



,(*) 



E {j) Ef E {k) Ef ] 



■ / dtie hK f ' t x 1 ' 1 S (Pi,p/) 

J 



(38) 



where = {Ej /c, pj) and similarly for p^\ Moreover: 

1 /■+# f+% e -(x!-x 2 ) 2 /44 



(39) 



(now we have explicitly indicated the integration volume). In order to compute S(pj,py), we 
change integration variables as follows: 



y = (xi+x 2 ) and x = (x : - x 2 ) , 



(40) 



and use the relation: 

H-f /■+# j_ /•+! /-+(X--*) 

/ etoi / dx 2 f (x 1 ,x 2 ) = - dx dy[f(x,y) + f(-x,y)] 

2 JO J-(L-x) 

Accordingly, we have: 



L „-x 2 /4r^. 



2 cos 



- ( P/ - Pi ) • X 



i=i 



Let us now take the limit L — > oo, which amounts to making the replacement: 

+0O 



1 



Yl ~^ I dp f alld J3 

P/=-°° 



1 



3 ' 



(41) 



(42) 



(43) 



In this limit, the term X{/L gives a vanishingly small contribution. Therefore we can write: 

,2 1a -2 



S (Pi,P/) 



1 



e£x- 



1 



(P/-Pi) 2 »-C 



(2vr/i) 3 7-oo (V^rcf i 2 ^ h ? ' 

The time integral in Eq. (I38p is trivial, and one arrives easily at the formula: 



(44) 



irijirikC 



u {p { f ] \ s f) u (Pi" ,Si)u(p\ K> ,Si)u(p [ f> \s f ) 



S3) 



(k) 



Ik) 



Im k 



E U) E U) E (k) E (k) 



(45) 
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As it is shown in the Appendix B, the integrating function (except for the Gaussian term) changes 
slowly within the region where the Gaussian term is appreciably different from zero. Therefore 
we can approximate it with the value it takes in the center of the Gaussian (where P/ = p.) and 
bring it out of the integral. Taking into account that u(p, s.) u (p, Sf) = S SiSf and performing the 
integration of the Gaussian part, Eq. (I45p takes the very simple expression: 



i m (vs-t)= rr^ m i mkC (46) 



(2) _____ 

We now turn our attention to the term I- (p_,s.;t) in Eq. ([37]) . Substituting Eq. (j33|) in 
Eq. ([37]) . we have: 

/ (2 V 8-t)--<y ^ V H^iJ- [dx, fdzz ^^^ c-^-^-^ (47) 

The spatial integrals are equal to 5(p, p.) defined in Eq. ([39]) . Performing the same type of 
calculation as before, and taking the limit L — > oo, one arrives at the result: 



/> <<P„ «!)- -7^/ rfp^-^e"^. (48) 

Once again, one can show that the integrating function (Gaussian term apart) varies slowly within 
the region where the Gaussian function is appreciably different from zero. Therefore one can bring 
this function out of the integral, fixing its value at the center of the Gaussian (p = p_), and perform 
the Gaussian integration. The final expression is: 

o\ 7m ■ rn^c 4 t 7T 3//2 

if (Pi,si;t) = "^4(7)7^3 V- (49) 
J 2 E W (2vr) d r 6 c 

Having computed explicitly all the terms of Eq. (136j) , we can turn our attention to Eq. (I34p . It 
is convenient to split the sum of Eq. (|34|) in one part with k = j and the another part with k ^ j: 



fc=l 
n 

+ _C [VakUpkV a jV} j P kj ( Pi , a.; i) + U^U^U^U^P,-* (p., _.;*)] . (50) 



Using the symmetry relation jjj^ = i^*, which implies that P&j (p., Sj, - 1) = P* k (pj, s.; i), one can 



rewrite Eq. (|50|) as follows : 

n n 

P Q ^ = ^ U* fc U^ fc U Qfc U^ + J2 2Re [U* fc U^ fc U Qi U^.P fcj (p., S< ; t)] , (51) 



fc=l fc=2 
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where we have exploited the identity: P^k (pi,sf,t) = 1. In the physically interesting case where 
the mixing elements U a k are real, Eq. (]5ip takes a very simple expression: 



P 



fe=l 



UafcU^fcUafcU^fc + 2. UafcU^fcUajU^j [1 — £jjfct] COS 



with: 



6fe 



167r 3 / 2 r^ 



v 7 ^ 



7 



2 4 

rrijC 



2^4 



mfe i?( fc ) 



167r 3 / 2 r^,mgC 4 



a) 



(fe) 



(52) 



(53) 



As a check, one can easily see that the probability is conserved, i.e.: 

1. 



(54) 



This is the result we wanted to arrive at. It shows that, also in collapse models, the number of 
particles is conserved, but the oscillations are damped according to Eq. (i52j) . by a factor equal 
to [1 — £jki\. This is in perfect agreemento with well established results concerning the effect 



of decoherence on oscillatory systems like those here considered |33j ]. Our calculation gives an 
analytical expression for the damping rate as predicted by the mass-proportional CSL model. 
Note that the calculation has been carried out to second perturbative order, which means £jkt <C 1. 
Therefore, the fact that the probability in Eq. (|52p becomes negative for ^j^t > 1 is of no concern, 
because this range of times goes beyond the limits of validity of the present result. Actually, one 
can try to stretch the above result beyond the second perturbative order, and guess the following 
expression for the transition probability: 



(El 



(fc) 



(55) 



= ^2 UafcU^fcUafcUgfc + ^ U^U^U^-U^ e ^ kt COS 
k=l k^j 

guments, which fix the form that the damping terms coming 



This is in agreement with general ar 
from decoherence effect should take 



37 



39]. 



The above results can be easily generalized to oscillatory systems, which decay in time. On 
the phenomenological level, one takes the decay into account by adding an imaginary term to the 
Hamiltonian: 



H 



h — r. 

2 



(56) 



The calculation remains unaltered, and one arrives at the final result: 

r (k) 



k=l 
n 



+ ^ UafcU/3fcUajTJ/3j ^ _ g - 



r( fc )+r(j) 



2 cos 



fc=2 



E^)t 



(57) 



In general, spontaneous collapses and decoherence are described by similar master equations. 
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cosmogenic 


solar 


laboratory 


E(eV) 


10 19 


10 6 


10 10 


too 


3.15 x 10 18 


5 x 10 2 


2,13 x 1(T 2 




2.31 x 1CT 55 


3.66 x 10~ 45 


1.56 x 1(T 57 



TABLE I. For our analysis we considered three types of neutrinos: cosmogenic, solar and laboratory neutri- 
nos. For each type, the table shows: the typical order of magnitude of the energies (first line), the time of 
flight (second line) and the damping factor as predicted by the mass proportional CSL model (third line). 

which generalizes Eq. ([52]) to decaying particles. 



V. THE CSL PREDICTION FOR NEUTRINO OSCILLATION 

Having derived the general formula, we now give a quantitative estimate the CSL effect on 
neutrino oscillations, by computing the damping factor in Eq. (I55p . We consider the stronger value 
7 ~ 10 _22 cm 3 s _1 for the collapse parameter, suggested in [15(. Substituting the numeric values 
of the constants in Eq. (|53|) . one finds that: 

djt ~ (7.33 x KT^s^eV 2 )-^ . (58) 

Here we have taken the largest possible squared mass difference m\c A — m\c A = 7.59 x 10~ 5 eV 2 j^j, 
where m\ and m-i are respectively the first and the second mass eigenstate and we have considered 



the ultra-relativistic approximation = \/pfc 2 + m^d 1 ~ p^c = E, and the same approximation 



(k) 

for Ef J . The energy E and the time t are free and depend on the nature of the neutrinos under 
study. Neutrinos detected in laboratories have mainly three origins: cosmogenic, solar and those 
produced in labs. Table I displays the typical values for the energy and time of flight (for simplicity, 
we assume that neutrinos travel at the speed of light) for these three types of neutrinos. The 
magnitude of the damping of the oscillations, as predicted by the mass proportional CSL model, 
has been evaluated using Eq. (|58j) . As we can see, in all three cases the CSL damping effect on 
neutrino oscillations is very small, and undetectable. The main reason is that the masses here 
involved — those of neutrinos — are very small, thus hampering the collapse mechanism. 
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VI. COMPARISON WITH THE DIOSI-PENROSE GRAVITY-INDUCED COLLAPSE 

MODEL 



As mentioned in the introduction, the damping of neutrino oscillations due to gravitational 
collapse, as described by the Diosi- Penrose model, was first studied in [23|]. In this work, the 
author argued that the decaying factor (the analog of for the CSL model) has the following 
form: 



G — 



D 



AEif (I/) dL' 



(59) 



where D is the distance such that AEq* (D) = 0, L is the distance traveled by the neutrino, and 

I r Mb#M^w' , (60) 



AE^ k 



where £ is a parameter that we will set equal to —G, with G the gravitational constant, like in the 
original paper by Penrose 241 ] . Moreover, pi (r) and p2 (r) are the two mass distributions, one for 



each different mass eigenstate. Since these two distributions travel at different velocities, AEt' k 



has a dependence on the traveled distance V (see Eq. (14) of [23J). Following the computation 



done in 



G 



23] and keeping all constants explicit, one finds: 



Ap? ~ 87T — 

G he 



3 (rrij + m-fc) ft? ra^ra^E ( 6 (raj + m^) irft?c 



L 



, „ (61) 
5Gf 2nftc \ hrrijmkGFE J 

where Gf is the Fermi constant and mj,rrik are the neutrino masses. In the case of cosmogenic 
neutrinos which have an energy of about E = 10 19 eV and travel a distance L ~ 10 25 m, the 
magnitude of the damping factor A^, fc lies between 1 and 10~ 2 , depending on the mass of the 
neutrino. Therefore, for the Diosi-Penrose model damping effect is by far stronger than that 
computed with the CSL model. This does not come as a surprise, because — as we stated in the 
introduction — it is well known that the model predicts, for single elementary constituents, a too- 



strong collapse of the wave function, which is incompatible with known experimental data 



2J. The 



reason is the following. It is clear that the model gives rise to divergences in the point-like limit, 
therefore one has to introduce a cutoff [3]. One way of doing it, is to consider elementary particles 
as spherical mass-distributions with a finite radius R. In [26J it was proposed to take R ~ 10~ 15 m, 
i.e. the nuclear size. However in j3] it was shown that also in this case the model is consistent 
with known facts (the energy increase of isolated systems, due to the collapse, is too large), and 
proposed a much larger radius, namely R ~ 10~ 5 cm, in order to restore compatibility. On the 



contrary, in [23J] the radius (aj, according to the paper's notation) R ~ Gprn/ft? ~ 10" 



-31±1 



m 



(where the uncertainty depends on the chosen value for the neutrino's mass) was considered. This 
cutoff is too small, therefore the result cannot be trusted. 
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VII. DECOHERENCE EFFECTS 

While traveling through the Universe and in particular through the atmosphere, neutrinos inter- 
act with the surrounding environment and scatter with other particles, mainly protons, electrons 
and other neutrinos. These interactions give rise to decoherence effects which also modify neutrino 
oscillations. Since protons interact with neutrinos only via the neutral weak current, both neutrino 
familieay are affected in the same way by this kind of interaction and neutrino oscillations are not 
modified. Unlike protons, electrons and neutrinos interact with the incoming neutrino both via 
neutral and charged weak currents: since these scatterings have different cross sections depending 



42] 



on the neutrino flavor, they contribute to decoherence 

A natural phenomenological estimate for the order of the decoherence rate is: A DEC ~ nva 
with v the relative velocity of the incoming neutrinos, n the density of the environmental leptons, 



43 



and a the relevant scattering cross section whose values are known in the literature 

<r„ e , e n 7 x KT^E^/GeVjcm 2 , (62) 

a Vlt>e ~ 10" 42 (E ! ,/GeV)cm 2 , (63) 

a VBtVe ~ 2, 8 x 10- 47 (E ! ,/GeV)cm 2 , (64) 

<Tu e ^ - 4 X 10- 48 (E ! ,/GeV)cm 2 . (65) 

The average density of electrons in the outer space and in the atmosphere are respectively n° UT ~ 
1/m 3 and n^™ ~ 2 x 10 26 /m 3 , while the average density of neutrinos is about n u ~ 10 8 /m 3 
everywhere. Assuming v to be the velocity of light in vacuum, then one finds: 

A Sec ~ KT 43 (£/eV) Hz , A£™ „ 10 -20 (E/eV) Hz. (66) 

with Aqec^™ 1 the decoherence rate in the out-space [atmosphere]. Neutrinos travel through the 
atmosphere within ~ 1CU 4 s. Using this data with the time-of-flights and energies in Table I, for the 
decoherence damping factor of the cosmogenic neutrinos, one finds: ~ 10~ 5 . For solar neutrinos 
instead, one gets: ~ 10 -18 , thus the damping of solar neutrino oscillations is hardly detectable, in 



agreement with experimental results 



4a 



46|]. 



This estimate shows that, since environmental decoherence on neutrino oscillations is much 
stronger than the CSL collapse effect, these spontaneous collapse effects cannot be observed exper- 
imentally, even if the technology were sophisticated enough to reach such sensitivities. They would 
anyhow be masked by decoherence effects. Moreover, decoherence effects are not far away from the 



In the following, we will consider decoherence effects only on electronic and muonic neutrinos. 
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collapse effect predicted in 



23 ]. which — as we argued — is anyhow overestimated. Therefore, also 



gravity-induced collapse effects cannot be detected via neutrino oscillations. 
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APPENDIX A: APPROXIMATION IN THE CALCULATION OF lf ] 

In this appendix we justify the approximation we used in order to derive Eq. (|49p from Eq. (|48p . 
This amounts to proving that: 

1 f 1 (P-P.)"^ 1 n 3 ' 2 

To see this, we can rewrite the integral in polar coordinates and integrate over the angular variables: 

1 (p-Pi) 2 ^ 2vr h 2 f +oc , p iv-vif-h , 



h3 J E { p 3) " & J.^ I p 2 c 2 + m 2 c 4 



Let us introduce the dimensionless variable s — 



(P-Pi)rc . 



1 fdV^-e-'-^ = - JL- r ds ^ (69) 



lipiC S> hc/rc ~ 10 eV (the typical range of momenta of neutrinos is between 10 3 eV and 10 19 eV) 
we can disregard s both in the numerator and denominator, obtaining: 



If- 1 _ (p-p,)H vr /■+- _ s2 tt 3 / 2 

-tt / dp— jrre & ~ / (Ise = pr, (70) 

tfJ E^ t%ePJ-oo r%E[ 3) 



which is the desired result. 
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APPENDIX B: APPROXIMATION IN THE CALCULATION OF 1$ 



Here we justify the approximation we used to pass from Eq. (|45p to Eq. (|46p . We start with 
Eq. (|45]): 



s f J J E\ 3] 'Ef 'Ef' 'EJ> 



1(1*0") _ T?U) _ Rl(fc) _l F ( fc )^ 



. (jB w_ jB o)_ jB w +jB w )(2 ^ r 

where: 



(71) 



Ef = Jp)c 2 + my , u(p,s)= fY<t m / 2^ s ) and V7-7 = ^— (72) 
' J ^J2mc 2 (E p + mc 2 ) m 

The first part of the integrand: 

^uipf^fMpf^^lsMpf^Sf) (73) 



jr<J) P Ci) P (*0 
^/ ^/ s / 



is a composition of polynomial functions of pj , and we can safely assume that it does not change too 
much, within the range where the Gaussian function is appreciably different from zero. Therefore 
we can then take P/ = Pi; by using also the relation: 

u{p?\sf)u(pf\ Si ) = 5 Sf:Si , (74) 



Eq. (fTTj) becomes: 



(i) nijinkC 1 / , 1 — e h > 1 ' * ip/ pj r c 



XX th * J X j % r 

V * ' 

=1 

Now we have to focus our attention on the integral /, which contains an oscillating term that needs 
special care. As before, we write the integral in polar coordinates and perform the integration over 
the angular variables; moreover, we introduce once again the a-dimensional variable s = ^ Pf * ■ 



We have: 

h 3 r+oc , n \ e ig(s) _l 2 

. — w t ds[ — s + pi) 

where we have defined: 



/ = ^rT3-M ds ( — S + Pi ) 7T73 e ~ s > ( 76 ) 



<,(*) = i(£7< fc > -EW-Ef + E®)t 



ct 

rc 



(s + y) + a k - J(s + y) + aj - y/y 2 + a k + Jy 2 + a,- , (77) 
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with 

2 

a-j = T -^m)<? = (l(T 2 eV- 2 ) m 2 c 4 e y = j- Pi = (lO^eV -1 ) Pi c (78) 

Our goal is to show that g (s) does not vary appreciably, within the range where the Gaussian 
term is significantly different from zero, and can be approximated with g (0) = 0; in this way, the 
integral can be computed exactly. This kind of approximation is not obvious because the factor 
ct/rc in front of Eq. (|77p can be very big. 

In the ultra-relativistic limit, we approximate the particle's velocity with the speed of light. 
In this limit, y 2 3> a,-,a/c, and so we can expand the square roots in g(s) using the Taylor series 
\Jx + e = y/x + 2^?j and we get: 

. . ct . , s ct (a,- — ah) , 
g (s ~ — ( aj -a k ) =—¥ rr- 79 



In order for g(s) to remain small within the interval where the Gaussian term of Eq. (176j) is 
appreciably different from zero, we need: 

y + — ) > — (oj-Ofc) (80) 
s J r c 

In all physically interesting situations, the term on the right hand side of Eq. (|80p is bigger than 1; 
moreover s is of the order of unity, because of the Gaussian in Eq. (1 76 p . Inequality ()80p is verified 
if the following condition is true: 



V > \ —{aj-a k ). (81) 
V rc 

Typically, cosmogenic neutrinos have energies bigger than 10 18 eV and travel distances of at most 
10 9 light-years 23]. This means that, in the worst case, ct/rc ~ 10 32 while (aj — ah) = (10 -2 eV~ 2 ) 
(m 2 c 4 - mf.c 4 ) ~ (10" 2 eV _2 )(2 x 10" 3 eV 2 ) = 1(T 5 . So we must have y > 10 14 , that means 
p iC = y/(10 _1 eV _1 ) » 10 15 eV, which is satisfied. 

For atmospheric neutrinos, ct/rc is in the range 10 11 - 10 14 while the range of energies is 



between 10 _1 GeV and 10 4 GeV 



4l(. This means that, even in the worse case, the condition to 



check is y ^> 10 5 , which means p^c = y/iW^eV' 1 ) > 10 6 eV. This is also satisfied. 
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